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' Abstract. This paper addresses the question of why metal-deficient barium stars are not yellow symbiotic stars 

0^ ■ (YSyS). Samples of (suspected) metal-deficient barium (mdBa) stars and YSyS have been collected from the 

1/^ ■ literature, and their properties reviewed. It appears in particular that the barium nature of the suspected mdBa 

' stars needs to be ascertained by detailed abundance analyses. Abundances are therefore derived for two of them, 

O ! HD 139409 and HD 148897, which reveal that HD 148897 should not be considered a barium star. HD 139409 is 

a mild barium star, with overabundances observed only for elements belonging to the first s-process peak (Y and 
I Zr). It is only moderately metal-poor ([Fe/H] = —0.4). The evidence for binarity among metal-deficient barium 

. stars is then reviewed, using three different methods: (i) radial- velocity variations (from CORAVEL observations) , 

CL(' (ii) Hipparcos astrometric data, and (iii) a method based on the comparison between the Hipparcos and Tycho-2 

JL I proper motions. An orbit is obtained for HIP 55852, whereas evidence for the (so far unknown) binary nature 

2 ■ of HIP 34795, HIP 76605, HIP 97874 and HIP 107478 is presented. No conclusion regarding the binary nature 

' of HIP 11595, HIP 25161 could be reached. Two stars with no evidence for binarity whatsoever (HIP 58596 and 

, BD -1-3° 2688) are candidates low-metallicity thermally-pulsing asymptotic giant branch stars, as inferred from 

IL" ' their large luminosities. The reason why mdBa stars are not YSyS is suggested to lie in their different orbital period 

distributions: mdBa stars have on average longer orbital periods than YSyS, and hence their companion accretes 
matter at a lower rate, for a given mass loss rate of the giant star. The definite validation of this explanation 
' 5^ ' should nevertheless await the determination of the orbital periods for the many mdBa stars still lacking periods, 

in order to make the comparison more significant. 
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1. The problem 

Our understanding of the link between chemically-peculiar 
red giants like barium stars or CH stars, and yellow sym- 
biotic stars (YSyS) has made substantial progre ss in the 
last decade (see the reviews 

by|jorisse3liQQ333), 

mainly 

with the realisation that likely all yellow symbiotics (i.e., 
involving a giant of spectral typ e G or K as primary com- 
poncnt) involve barium stars ijSmith et all Il996l Il997t 
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IPereira fc Porto de Mellolll997HPereira et al.l . ll998() . The 

metal-deficient nature of the giant is a key factor, because 
it implies a rather large luminosity for the giant. Since 
evolutionary tracks of metal-deficient stars are shifted to- 
wards the blue, metal-deficient giants of spectral type 
K must lie on the uppe r part of the (asym ptotic) giant 
branch (see Fig. 11 of ISmith et all Il996|) . where they 
suffer strong mass loss. If such metal-deficient K giants 
are in binary systems, their strong wind will interact with 
the companion and trigger symbiotic activity. 

The other facet of this problem, namely whether all 
metal-deficient barium stars are symbiotic stars, is not 
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yet fully answered. The present paper offers a first step in 
that direction. 

We have collected a list of candidate metal-deficient 
barium stars and have assembled new observations to 
check (i) whether these stars are binaries, (ii) whether 
they are barium stars and (iii) whether they exhibit sym- 
biotic activity. 

2. The samples 

Before discussing metal-deficient barium stars, it is useful 
to first summarize the properties of YSyS, to which metal- 
deficient barium stars may be compared. 

2. 1 . Yellow symbiotic stars 

All known YSyS are listed in Table which shows that 
all the stars studied so far exhibit the barium syndrome. 
YSyS with a stellar infrared continuum (s-typc, as op- 
posed to the dusty d'-type; see below) are clearly metal- 
deficient objects, as revealed by their low metallicities and 
high space velocities (CD —43° 14304 may be an exception; 
however, it is of spectral type K7, and should perhaps not 
be included in the family of YSyS). The presence of the 
barium syndrome among a family of binary metal-deficient 
stars fully supports the commonly accepted hypothesis 
that the s-process is more effi cient at low metallicities 
(|ClavtonL Il988t I Jorissenl l2003a|) . s-Typc YSyS, with their 
metallicities lower than classical barium stars, may be ex- 
pected to be, on average, rnore l uminous than the latter 
(see Fig. 11 of lSmith et al.[ll99(T[ comparing the luminos- 
ity function of Pop. I and Pop. II K giants). This is a direct 
consequence of the fact that evolutionary tracks shift to- 
wards the blue in the Hertzsprung-Russell (HR) diagram 
as metallicity decreases, as shown in Fig.^. Fig.^ con- 
firms that the YSyS AG Dra and BD -21°3873 are indeed 
more luminous than classical barium stars. This difference 
in the average luminosity - and hence mass-loss rate - of 
the two populations thus explains why YSyS, despite host- 
ing a K giant, exhibit symbiotic activity whereas barium 
stars do not. The larger mass-loss rates for the cool com- 
ponents of s-typc YSyS - as compared to Ba stars - may 
be inferred from the comparison of their IRAS [12] — [25] 
color indices, which reflect the amount of dust present in 
the system: ([12] — [25])Ra. < 0.1, as compared to 0.45 
for AG Dra JSmitdi et allll99(i^ . lMiirset et all ((l99l and 
iDrake et al.l (|l9871) provide direct measurements (or up- 
per limits) for the mass loss rates of AG Dra and of Ba 
stars, respectively, which confirm the above conclu sion. 

Y SyS with a dusty infrared contin uum (d'-type; lAUenl . 

Il982l: ISchmid fc NiissbaiimeJ. Il993l) differ from their s- 
type counterparts in several respects (Table ^ : they host 
a complex circumstellar environment (including cool dust, 
bipolar outflows, extended optical nebulae or emission-line 
spectra closely resembling those of planetary nebulae) , the 
cool components have early spectral types (F to early K), 
they are often fast ro tators (with the possible exception 
of M 1-2 =V471 Per; iGranier B(^ll98lh and, finally. 



they belong to the galactic disk unlike s-type YSyS which 
belong to the halo. 

All these arguments suggest that the hot component 
in d'-type SyS has just evolved from the AGB to the 
WD stage. The rather cool dust (|Schmid fc Nussbaumeil 
Il993t) is a relic from the mass lost by the AGB star. The 
optical nebulae observed in d'-type SyS arc most likely 
genuine planetary nebulae rather than the nebulae as- 
sociated with the io nized wind of the cool component 
l|Gorradi et all. Il999t) . This is esp ecially clear for AS 201 
which actually hosts two nebulae 1 Schwar^ll99lh : a large 
fossil planetary nebula detected by direct imaging, and a 
small nebula formed in the wind of the current cool com- 
ponent. Finally, the rapid rotation of the cool component 
has likely been caused by spin accreti on from the former 
AGB wind like in WIR RING systems ()jeffries fc Stevensl 
.1996: i.Torisscnl l2003lJ ). The fact that the cool star has 
not yet been slowed down by magnetic braking is another 
indicatio n that the mass transfer occurred fairl y recently 
(Theuns et al.Ul996|) . ICorradi fc Schwar j l|l997|) obtained 
4000 y for the age of the nebula around AS 201, and 
40000 y for V417 Gen. 



2.2. Metal-deficient barium stars 

Metal-deficient barium stars (with metallicities in the 
range —1.1 to —1.8, comparable t o that of YSyS) 
were i dentif ied bvlLuck fc Bqndl l)l99l|) , iMennessier et alJ 
((l997t) and IZacs et alJ l|200fl^ ~and occupy the same re- 
gion of the HR diagram as YSyS (Fig.^). The question 
thus arises why metal-deficient barium stars are not SyS. 
Different answers must be sought, depending upon their 
absolute visual magnitudes My. The most luminous sys- 
tems, with Mv < —2, are likely located on the thermally- 
pulsing AGB, so that their Ba syndrome may be explained 
by internal nucleosynthesis. They thus should not be bi- 
naries, and therefore cannot be SyS! HD 104340 (open 
circle in Fig. Hb), a metal-deficient Ba star studied by 
lJunaueira fc Pereiral l|200l|) . and BD -f03°2688 (Table 
provide good illustrations of this situation, since they both 
lie above the TP-AGB threshold and CORAVEL radial- 
velocity measurements spanning several years do not re- 
veal any clear orbital motion (Figs. |5| and ^1 as well as 
Scct.EJ. 

The less luminous and warmest among metal-deficient 
Ba stars, clumping around Afy ~ +1 in the HR dia- 
gram, are also sometimes classified as CH stars (crosses 
in Fig. ^). They are not losing mass at a large enough 
rate to trigger any symbiotic activity, as revealed by thei r 
small [12] - [25] color indices (< 0.3: ISmith et al.lil996|) . 

Finally, at intermediate luminosities (—2 < My < +1) 
where YSyS are located, metal-deficient Ba stars are not 
luminous enough to be TP-AGB (hence they should be bi- 
naries), but yet their mass loss rates must be large enough 
to trigger symbiotic activity (provided that the orbital 
separation is not too large, since it is the mass accretion 
rate by the compact companion which is in fact the key pa- 
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Table 1. The barium syndrome among YSyS. The spectral type of the cool component is taken from Miirset & Schmid 
(1999), or references therein. In column labelled 'nebula', 'y' means that an optical nebula has been detected, and 'PN' 
that, based on its emission line spectrum, the star has traditionally been included in planetary n ebulae catalogues, 
even t hough no optical nebula may be visible. The column labelled 'P' lists the orbital period, from lMiirset fc SchmidI 
lll999t) 



Name 


Sp. Typ. 


[Fe/H] 
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[Ba/Fe] 


V sin i 


nebula 


P 


Ref. 








(km/s) 


(°) 




(km/s) 




(d) 




d'-type 


V417 Cen 


G8-K2 


~ 0.0 




-1 


0.5 


70 


y 


247 


(5,11) 


HDE 330036 


G5 


0.02 


-14 


+4 


0.88 


100 


PN 




(14,17) 


=Cn 1-1 




















AS 201 


G5 


0.07 




+7 


0.63 


25 


y 




(12,17) 


V471/V741 Per 


G5 


? 


-12 


-9 


> 




PN 




(2) 


=M 1-2 




















St He 190 


G5 


0.0 


~ 10 


-35 


~ 0.5 


100 


bip. outf. 




(10,13) 


Wray 157 


G5 


? 
















Hen 1591 


< K4 


? 
















s-type 


UKS Ce-1 


C4,5Jch 


? 




+2Q 


> 








(6) 


S 32 


C1,1CH 


? 


-h325 


-30 


> 






612 


(6,14) 


Hen 2-467 


KO 


-1.1 


-109 


-12 


+0.8 




n 


478 


(4,16) 


BD-21°3873 


K2 


-1.1 


+2M 


+37 


+0.5 




n 


282 


(3,15,16) 






-1.3 






+0.3 








(9) 


AG Dra 


K2 


-1.3 


-148 


+41 


+0.5 




n 


554 


(8,16) 


CD -43° 14304 


K7 


-1.4 


+27 


-41 


? 






1448 


(7,18) 



-1.0 



Chemical evolution of the Galaxy 
< 0.2 



(1) 



References: (1) Edvardsson et al., 1993, A&A, 275, 101 (2) Grauer & Bond, 1981, PASP, 93, 630 (3) Pereira et al., 1997, AJ, 
114, 2128 (4) Pereira et al., 1998, AJ, 116, 1977 (5) Pereira et al., 2003, In: Symbiotic stars probing stellar evolution, R.L.M. 
Corradi, J. Mikolajewska, T.J. Mahoney (eds.), Astron. Soc. Pacific Conf. Ser. (San Francisco), p. 85 (6) Schmid, 1994, A&A, 
284, 156 (7) Schmid et al., 1998, A&A, 329, 986 (8) Smith et al., 1996, A&A, 315, 179 (9) Smith et al., 1997, A&A, 324, 97 (10) 
Smith et al., 2001, ApJ, 556, L55 (11) Van Winckel et al., 1994, A&A, 285, 241 (12) Schwarz, 1991, A&A, 243, 469 (13) Munari 
et al., 2001, A&A, 369, LI (14) Schmid & Nussbaumer, 1993, A&A, 268, 159 (15) Munari & Patat, 1993, A&A, 277, 195 (16) 
Corradi et al., 1999, A&A, 343, 841 (17) Pereira et al., 2005, A&A 429, 993 (18) The metallicity is from Pereira, priv. comm. 



rameter; see Sect. Eland Jorissen 2003a). It is thus of great 
interest to check (i) the Ba nature of those metal-dcficicnt 
stars with intermediate luminosities, (ii) their binary na- 
ture, and (iii) their suspected symbiotic activity. The first 
two issues are addressed in Sects. 13 and 0] respectively. 

As far as a possible symbiotic activity is concerned, 
there is no indication from their photometric U — B and 
B — V indices that the mctal-dcficient stars in Table |21 
have a strong blue continuum which could betray their 
symbiotic nature. It is thus very likely that none among 
these stars is a full-fledged symbiotic star. No signature 
of weak symbiotic activity (of the kind exhibited by some 



binary S stars; sec Fig. 1141 and IVan Eck fc Jorissenil200^ 
is observed in the Hq, line profile either (Fig. [21). 

3. Abundances 

The classification of the stars in Table |21 as metal- 
deficient Ba stars is subject to caution, as it does not 
rely on spectroscopic abundance analyses, but rather on a 
maximum-likelihood assignment based on kinemat i c, spa- 
tial and luminosity properties |Mennessier et all . Il997|) 
for barium stars from the list of Nevertheless, 
when a metallicity determination is available, it confirms 
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Fig. 1. Left panel ( a): Evolutionary tracks oflSchaller et alJ l)l992|) compared with the locations of classical barium stars 
(stars labelled 'G' in lMennessier et all <|l997f) : open dots) and the yellow SyS AG Dra (filled circle), BD -21°3873 (filled 
triangle) and Hen 2-467 (filled square). The bolometric magnitudes were take n from the re f erence s listed in Tabled 
Thes e bolometric magnitu des were combined with bolo metric corrections from lBessell et all l)l998|) and B — V indices 
from iMunari et al ] l|l992l) and lMunari fc BusonI l)l992|) to yield the absolute visual magnitudes. 

Right panel (b): Same as (a) but for YSyS (filled symbols as in the left panel) and metal-deficient b arium stars [open 
triang les: stars flagged as 'H' bv iMennessier et al. l)l997|) : crosses: CH stars also flagged as 'H' by iMennessier et alJ 



small open squares: additi onal metal-deficient, s-proc ess-rich stars from Table |21 open circle: HD 104340, 
large open square: HD 206983 from ljunou eira fc Pereiral ll200lM . The dashed horizontal line represents the luminosity 
(-Wboi = —3, corresponding to My ~ —2) at the first thermal pulse in a 1 Mq AGB star of metallicity [Fe/H] =—1.8 
according to iLattanziol l)l99l[l . 



the metal-deficient nature of the object (see Table EJ. 
HD 139409 is an exception, though, since detailed spectral 
analyses reveal that it is neither metal-poor nor strongly 
enriched in s-process elements (see below). It may never- 
th eless be hoped that the m etal-deficient assignment made 
bv IMennessier et al.l 1 1997t) is valid in all the other cases. 
Re gard i ng th e Ba nature of these stars, it is known that 
thelLiil l(l99lh catalogue of bari um stars, from wh i ch the 
sample of barium stars used by IMennessier et al ] ^ll997^ 
was drawn, is contaminated by many non-b arium stars 
l)Griffin fc KeenanL Il992t Ijorissen et all Il996|) , especially 
among those stars having a Ba index smaller than 1. 



It would therefore in principle be necessary to re- 
evaluate the Ba nature of all the stars listed in TableO So 
far, spectra could be obtained for two of them, HD 139409 
and HD 148897, which are discussed in detail in the 
present section. Of these, only HD 139409 appears to be 
a mild barium star, thus confirming the suspicion about 
theiLii ((199lil catalogue expressed above. 



3.1. The case of HD 148897 and HD 139409 

HD 148897 (= HR 615 2) has been tagged as a 'hkely 
marginal barium star' b vlBovle fc McClurg 1 197lj) and as 
a marginal CH star by IVilhu et all l|l977t) . It therefore 
found its way i nto th e barium-star catalogs of iLii et alJ 
(|1983|) andO l|l99l|) . as well as the catalog of CH an d 
metal-deficient Ba stars of lSleivvte fc BartkeviciusI l)l990|) . 
The star was clas s ified as G8.5III CN-2 Fe-1 CH-1 by 
iKeenan fc McNe"!]! ljl989l) . and this classification as CN- 
and CH-weak contradicts the earlier assignments. 

A detailed abundance analysis has thus been per- 
formed to clarify the situation, and its results ar e com- 
par ed with prev ious analyses bv iKvrolainen et aP 
a,nd iLuckl lll 99lh in Sect. \rr^ 

HD 139409 (= HIP 766 05) has been class i fied 
a marginal bariu m star by iMacConnell et al) lll972l) . 
as G5 HI Bal by lYamashita. fc Norimotol l|l98l|) and as 
KOIII/H Ba 0.5 bv lLuNi 99lf> " 



3.1.1. Observations 

A high-resolution spectrum of HD 148897 was ob- 
tained using the Coude Matrix Echelle Spectrometer 



as 
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Table 2. Stars classified as metal-deficien t barium stars by Mennessier et al The absolute magnitude My is a 

maximum likelihood estimate obtained bv iMennessier etaL except for the additional stars where it is derived 

from a straight inversion of the Hipparcos parallax (BD +75°348); from a spectroscopic estimate of gravity and an 
educated guess for the mass (BD-|-3°2688) or from a fit to the M 92 isochrone (CS stars). The column labelled 'Ba' 
indicates whether detailed chemical analyses have confirmed the Ba nature of the star. 



HIP 


HD/DM 


Mv 


B-V 


[Fe/H] 


Ref 


Ba 


Ref 


Rem 


4347 


5424 


-0.6 


1.14 












11595 


15589 


-0.5 


1.15 


—0.7 


9 










— Zi ZZoo 


-i.i 


l.io 












29740 


43389 


-3.5 


1.48 












34795 


55496 


-1.7 


0.96 


-1.55 


5 


y 


5 




58596 


104340 


-1.9 


1.36 


-1.72 


3 


y 


3 




69834 


123396 


-0.9 


1.19 












76605 


139409 


-0.5 


1.04 


-0.42 


2 


mild 


2 


AIv = 1.5 is derived from spectroscopic logg 


80843 


148897 


-1.8 


1.27 


-1.0 

-0.62 

-1.16 


2 
7 

8 


n 


2 




97874 


187762 


-1.2 


1.10 












107478 


206983 


0.2 


1.64 


-1.43 


3 


y 


3 





Additional stars [from refs. (1), (5) and (10)] 



43042 +75°348 -0.4 : 1.17 -0.8 1 y 1 

55852 +4°2466 ? 0.60 -1.85 5,11 y 4,5 

-h3°2688 -5.0 1.22 -1.42 5,11 y 5 

CS 22942-019 0.74 -2.67 10 y 10 

CS 22948-027 -1 0.93 -2.60 10 y 10 



Ref erences to the Tab le: ^ ^ ^ ^ 

(1) IZacs et al.N2000l) ( 2) This work (31 lJunaueira fc Pereiral j200l|) (4) iBurris et alJ tOOd) (5) iLuck fc Bone 
IJorissen et aLrjl998tl f7l iKvrolainen et alJ lll98el) (8 ) lLuckl diggiT ) fgl Barbuv et all Jiggl flO) |Preston fc Sneder 
B-V from the Tvcho-2 catalog ijHag et all l2000b^ 



CMAESTRO: iMusaev et"ail Il999j) delivering a resolving 
power of 45 000 and installed on the 2-m Zeiss telescope 
of the Terskol Observatory (located in Northern Caucasus 
at an altitude of 3100 m). The spectrometer is equipped 
with a Wright Instruments CCD detector with 1242 x 
1152 pixels (22.5 x 22.5 /im). A total exposure of 1800 s 
was taken on February 18, 2003. The spectrum covers the 
range 365 to 1020 nm spread over 85 spectral orders. 

A high-resolution spectrum of HP 139409 was ob - 
tained on the HARPS spectrograph ijMavor et al.Ll2003() . 
delivering a resolution of 115 000 and installed on the 
ESO La Silla 3.6 m telescope. A total exposure of 200 s 
was obtained by Xavier Bonfils on May 31, 2004. In or- 
der to check equivalent widths delivered by HARPS, a 
spectrum of the standard star Arcturus was obtained as 
well by Fabicn Carrier. The HARPS spectra were re- 
duced by the observers using standard pipeline processing. 
Equivalent widths for the same set of lines as those studied 
in HD 139409 have been measured by one of us (L.Z.) in 
the HARPS spectrum of Arctur us and compared to those 
from the Arcturus spectral atlas ijHinkle et al.l 12000}) . The 
agreement between the two sets of equivalent widths is 
excellent (Fig. OJ, thus qualifying HARPS for abundance 
analyses. Spectra around the A614.172 nm Ba II line are 



shown in Fig.Blfor t he two analyzed stars and for Arcturus 
llHinkle et al.n200nl) . 

3.1.2. Reduction and analysis 

HD 148897 

The reduction of the CCD frames (subtraction of bias, 
dark and scattered light, flat flclding, extraction of cchcUc 
orders and wavelength calibration) was per formed with 
the DECH20T software l)Galazutdinovlll992|) . More than 
500 weak to medium-strong atomic lines, free of blends, 
were identified and their equivalent widths were measured 
with the DECH routines. The radial velocity was mea- 
sured using a large number of symmetric absorption lines. 
The observed velocity has been brought to the heliocentric 
system by adding +22.7 km s~^ . The mean heliocentric ra- 
dial velocity for HD 144897 was found to be +16.7 km s^^ . 
The atmospheric parameters for that star cannot be de- 
rived from photometry, since the standard temperature 
calibrations refer to stars of normal chemical composition. 
To obtain a colour-independent estimate of the tempera- 
ture, a spectroscopic temperature has been derived from 
the excitation equilibrium of Fe I lines. The surface grav- 
ity logg was determined using the Fel/Fell ionization bal- 
ance, whereas the microturbulent velocity Vt was derived 
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Fig. 2. Ha profiles for two candidate metal-deficient Ba 
stars (HD 139409 and HD 148897; see however the discus- 
sion in Sect. 131), compared to that of the K giant Arcturus. 
The wavelength scale has been corrected for the stellar ra- 
dial velocity. The profiles show no indication of symbiotic 
activity. 



by forcing the abundances of individual Fc I lines to be 
independent of the reduced equivalent width. The result- 
ing atmospheric parameters for HD 148897 are as follows: 
Teff = 4350 K, logg = 1.0 (cgs), and Vt = 2.0 km s'^. 
An independent determ ination of the surface gr avity, us- 
ing the above Tcff vahie. lMennessier et alJ (^23) absolute 
visual magnitude (Afy = —1.8), a bolometric correction of 
—0.5 and a mass of 1 Mq, yi elds log q = 1 -1, in agreement 
with the adopted value (see iLuckL Il99ll for a discussion 
of the discrepancy usually observed between the photo- 
metric and spectroscopic gravities). Comparison with at- 
mospheric parameters from the literature is presented in 
Table El 

The abundance analysis has been performed with 
the standard LTE line analysis program WIDTH9 de- 
velo ped by Kurucz. The mo del atmospheres were taken 
from lGustafsson et al ](|l97i). The synthetic spectra were 

! :enerated using the spectral synthesis code STARSP 
Tsvmbal Il996^ . Oscillat or strengths have be en taken 
from the VALD database ijPiskunov et al.l.ll995|) . The re- 
sulting ab undances (norma lized by the solar-system abun- 
dances of ICrevesse fc Sauvak 1199?!) are listed in Table 



from which it may be concluded that HD 14 8897 appears 
to be a rather typical mctal-dcficient star (jMcWilliaml . 
Il997l) . and should certainly not he considered as a (metal- 
deficient) barium star. 

Interestingly enough, there arc no indications whatso- 
ever from the Hipparcos and Tycho data (applying the 
methods described in Sect. 0} that this star is binary, in 
agreement with the fact that it is not a barium star. 

HD 139409 

The effective temperature of HD 139409 has been derived 
from the excitation equilibrium of Fei, Til and Cri lines 
(see Fig. for Fe) . The surface gravity log g was deter- 
mined from the Fei/Feii ionization balance, and the mi- 
croturbulent velocity by forcing the abundances of individ- 
ual Fei, Til and Cri lines to be independent of equivalent 
width (see Fig. for Fe) . Although spectroscopic grav- 
ity and temperature determinations in late-type, metal- 
deficient stars are probably affected by non-LTE effects, 
these effe cts remain small when [Fe/ H] > —1.0 (see, for 
example, lAllende Prieto et al.Ul99'9(l . The stellar param- 
eters of Arcturus ([Fe/H]= —0.6, logg — 1.3) derived by 
the spectroscopic method (applied on the HARPS spec- 
trum) are in good agreement indeed with those derived by 
other (non-spectroscopic) methods. 

The resulting atmospheric parameters for HD 139409 
are as follows: Toff = 5000 K, log.g = 2.8 (cgs), and £_t = 
2.0 km s^^. The spectroscopic gravity, combined with a 
mass of 1 Mq, leads to My ~ -1-1.5. Thus our calculations 
indicate that HD 139409 is less luminous than predicted 
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Table 3. Averaged absolute ab undances (in the scale wh ere loge(H) = 12) and normalized abundances [X/Fe] for 
HD 148897 (relative to the Sun; |g revesse fc Sauvalll99i) . Th e standard d eviations a and the number of lines used 
in the analysis (n) are also given. Two sets of abundances from lLuc^ l|l99 J) are given: one is based on a 'photometric 
logg' (derived from the absolute magnitude, the effective temperature and an estimate for the stellar mass) and the 
other on the 'spectroscopic logg' (Fe I/Fe II ionization balance). 

This work Kirolainen et al. (1986) Luck (1991) Luck (1991) 

Teff= 4350 K Teff= 4360 K T^fi= 4100 K T^a= 4100 K 



X 




Vt 


log 3 = 1.0 
= 2.0 km s"^ 






lo 

Vt = 


gg = 1.5 
2.6 km s~^ 


phot 

Vt = 


logg = 1.6 spec. \ogg = 0.1 
1.7 km s~^ Vt = 1.7 km s~^ 


loE e(X) 


a 


n 


L / J W 






[X/HIr 

o 


[X/Felfii 

L / J W 






,11 


[X/H]o 

o 


^i.ii 








O I 


8.22 


0.20 


4 


-0.61 


H-0.41 


-0.58 


+0.04 


-0.50 


+0.50 


-1. lU 


+0.06 


Na I 


5.10 


0.08 


6 


-1.23 




0.21 


-1.11 


-0.49 


-1.36 




0.36 


-i.zD 




0.10 


Mgl 


6.92 


0.09 


6 


-0.66 




fO.36 


-0.59 


+0.03 


-0.54 


+0.46 


-U.DZ 


+0.54 


Al I 


5.61 


0.23 


4 


-0.86 


-f0.16 


-0.70 


-0.08 


-1.28 




0.28 


1 on 

-i.ZU 




0.04 


Si I 


6.91 


0.11 


23 


-0.77 


-fO.25 


-0.56 


+0.06 


-0.24 


+0.76 


-U.oy 


+0.57 


K I 


4.17 




1 


-0.95 


-f0.07 


- 


- 


- 












Ca I 


5.65 


0.11 


21 


-0.71 


+0.31 


-0.54 


+0.08 


-1.09 




0.09 


-u.y / 


+0.19 


Ca II 


5.48 




1 


-0.88 




f0.09 


- 


- 


- 












Sc II 


2.08 


0.12 


13 


-1.09 




0.12 


-0.86 


-0.24 


-0.75 




0.40 


-i.o / 




0.21 


Ti I 


4.20 


0.07 


56 


-0.82 


4-0.20 


-0.48 


+0.14 


-1.20 




0.20 


1 1 Q 




0.03 


Ti II 


4.22 


0.14 


8 


-0.80 


+0.17 


- 


- 


-0.51 




0.16 


-i.lD 


0.00 


V I 


3.11 


0.09 


23 


-0.89 




f0.13 


-0.56 


+0.06 


-1.24 




0.24 


-i.z4 




0.08 


V II 


3.01 




1 


-0.99 




0.02 


- 


- 


- 












Cr I 


4.59 


0.12 


27 


-1.08 




0.06 


-0.82 


-0.20 


-1.19 




0.19 


-i.Uo 


+0.08 


Cr II 


4.85 


0.08 


6 


-0.82 


+0.15 


- 


- 


- 












Mn I 


4.24 


0.10 


11 


-1.15 




0.13 


-1.00 


-0.38 


- 












Fe I 


6.48 


0.07 


180 


-1.02 






-0.62 


- 


-1.00 






-1.16 






Fe II 


6.53 


0.07 


17 


-0.97 






- 


- 


-0.35 






-1.16 






Co I 


4.11 


0.12 


21 


-0.81 


+0.21 


-0.70 


-0.08 


-0.84 


+0.16 


-1.10 


+0.06 


Ni I 


5.15 


0.11 


54 


-1.10 




0.08 


-0.62 


0.00 


-0.84 


+0.16 


-1.12 


+0.04 


Cu I 


3.31 


0.17 


3 


-0.90 


+0.12 






-0.64 


+0.36 


-0.80 


+0.36 


Zn I 


3.78 




1 


-0.82 


+0.20 






-0.39 


+0.61 


-0.92 


+0.24 


Y II 


1.21 


0.12 


6 


-1.03 




0.06 


-0.73 


-0.11 


-0.62 




0.27 


-1.14 


+0.02 


Zr I 


















-1.43 




0.43 


-1.43 




0.27 


Zr II 


1.91 


0.13 


3 


-0.69 




fO.28 


-0.37 


+0.25 














Ba II 


0.97 


0.05 


3 


-1.16 




0.19 


-0.66 


-0.04 














La II 


0.24 


0.25 


3 


-0.93 


+0.04 


-0.52 


+0.10 














Ce II 


0.26 


0.16 


5 


-1.32 




0.35 


-0.65 


-0.03 


-0.68 




0.33 


-1.21 




0.05 


Pr II 


-0.20 


0.13 


2 


-0.91 


+0.06 


















Nd II 


0.56 


0.13 


8 


-0.94 


+0.03 






-0.61 




0.26 


-1.15 


+0.01 


Eu II 


0.17 




1 


-0.34 




fO.63 






-0.40 




0.05 


-1.01 


+0.15 



© 



bv iMennessier et al The derived iron abundance 

is [Fe/H] = -0.42 [adopting log(e(Fe)) = 7.50]. 

The basic conclusion from the set of abundances listed 
in Table H and displayed in Fig. □ is that HD 139409 
appears to be a mild barium star.^ The s-process over- 
abundances observed in HD 139409, although quite mod- 
erate, are not much smaller than those observed in a yellow 
symbiotic star like BD -21°3873, which exhibit overabun- 
dances of s-process elements in the range 0.3 - 0.8 dex 

^ Interestingly, iPinsonneault et al.l I^S^) quoting a private 
communication from E. Luck, remark that Luck "is complet- 
ing an abundance analysis of HD 139409, which will show that 
the light s-process elements do have (~ +1.5 dex) enhance- 
ments, while heavier s-process element enhancements are much 
smaller (~ +0.3 dex)". Luck's study seems to have never been 
published. 



(ISmitheE3[T993). With [Fe/H] = -0.4, HD 139409 is, 
however, not as metal-deficient as the other stars consid- 
ered in this paper. 

4. Statistics of binarity among metal-deficient Ba 
stars 

The list of confirmed or suspected metal-deficient barium 
stars remaining after the screening process based on abun- 
dances (as described in Sect. is given in Table [3 The 
possible binary nature of those stars may be assessed using 
three different methods: 

— checking for radial-velocity variations; 

— checking for astrometric orbital motion, directly from 
the Hipparcos astrometric data; 
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Hinkle et al . (2000) 

ff] 
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Fig. 4. The spectra of the target stars in the region of 
the Ba II line at 614.172 nm. Measured lines are marked. 
Also shown is the spectrum of the standard star Arcturus 



-0.2 - 




1 2 3 4 5 

Lower Excitation Potential (eV) 



).2 



„ o 



00°, 



)OCO 8™ 
» o» 



• CPA 



o o 
• 

8 o 

o o 
%o°o 



20 40 50 

Equivalent Width (mA) 

Fig. 6. The abundances derived from the Fe I (open cir- 
cles) and Fe II (filled circles) lines are displayed as a func- 
tion of the line equivalent widths. The absence of a trend 
in these data is used to derive the microturbulent velocity 
of HD 139409. 



0,5 



-0.5 



Na 
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Fig. 7. Abundance pattern in HD 139409 (abundances 
derived from ionized species are represented by black dots, 
and from neutral species by crosses). 



Fig. 5. The iron abundances derived from Fe I (open cir- 
cles) and Fe II (filled circles) lines are displayed as a func- 
tion of the excitation potential for the lower energy level 
of the line. The absence of a trend in these data is used 
to derive the spectroscopic temperature of HD 139409. 



— checking for astromctric orbital motion, indirectly 
from a comparison of the Hipparcos and Tycho-2 
proper motions. 



These methods are now described in turn. 
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Table 5. Binary properties of confirmed or suspected metal-deficient barium stars. The column labeled 'Ba' indicates 
whether detailed chemical analyses have confirmed the Ba nature of the star. The columns labeled X^Hip-pxyc ^^'^ 
Prob(x^ > X^Hip-MTyc) provide the (and its associated probability) involved in the comparison of Hipparcos and 
Tycho-2 proper motions (see text for details). The column 'Bin.' has been set to 'y' if the first kind risk of rejecting 
the null hypothesis that the Tycho and Hipparcos proper motions are equal is smaller than 10%. The column labeled 
IAD indicates whether the signature of an orbital motion is present in the Hipparcos Intermediate Astrometric Data, 
according to the various tests described in the text. The column labeled a{Vr) provides the radial-velocity standard 
deviation, At and TV correspond to the time span and number of observations, respectively. The columns labeled P 
and 'Ref list the orbital period (when available) and the reference for the radial velocity and/or orbital data. The 
column 'Binary' gives the final binarity diagnostics. 



HIP 


HD/DM 


Ba 


Ref 


zu (mas) 


Xobs 


/iHIP — MTyc 

Prob(x" > XobJ 


Bin. 


IAD 


a{Vr) 
(km s^^) 


At 
(d) 


N 


P 
(d) 


Ref 


Binary 


4347 


5424 






0.22 ± 1.42 


1.05 


0.59 


n 


n 


2.32 


3306 


13 


1881 


6 


y 


11595 


15589 






2.03 ± 1.21 


1.63 


0.44 


n 


n 












? 


25161 


-27°2233 






0.89 ± 1.35 


0.23 


0.89 


n 


n 












? 


29740 


43389 






-1.25 ± 1.00 


4.54 


0.10 


y 


y 


3.85 


3350 


24 


1689 


6 


y 


34795 


55496 


y 


5 


2.44 ± 1.04 


0.02 


0.99 


n 


y 


0.73 


5121 


24 




2 


y? 


58596 


104340 


y 


3 


-0.97 ± 1.09 


3.39 


0.18 


n 


n 


1.64 


2587 


16 




2 


n 


69834 


123396 






1.73 ±0.86 


1.98 


0.27 


n 


n 












? 


76605 


139409 


y 


2 


5.51 ± 1.14 


4.98 


0.08 


y 


n 


0.66 


1478 


2 




2 


y? 


97874 


187762 






2.07 ± 1.53 


9.00 


0.01 


y 


y 












y 


107478 


206983 


y 


3 


3.75 ±1.86 


2.94 


0.23 


n 


y 












y? 



Additional stars [from refs. (1), (5) and (10)] 



43042 


±75°348 


y 


1 


1.02 ± 1.32 


0.12 


0.94 


n 


y 4.64 


1436 


33 


1042 


11 


y 


55852 


±4° 2466 


y 


4,5 


0.96 ± 1.83 


0.43 


0.81 


n 


n 3.38 


6191 


41 


4592 


2 


y 




±3°2688 


y 


5 










0.45 


986 


7 




2 


n? 




CS 22942-019 


y 


10 










3.42 


3274 


15 


2800 


10 


y 




CS 22948-027 


y 


10 










1.87 


2560 


13 


505 


10 


y 



References to the Table: 



(1) IZacs et al.l J2000D (2) This work (3) IJunnueira fc Pereiral j200lD (4) iBurris et alJ (1200(1) (5) iLuck fc Bonl Jl99lD (6) 
lUdrv et all lll998h flO) [Preston fc Sneded i200ll) flD lZacs etaO tOO^ 



4.1. Radial-velocity variations 

Several stars from Table (namely HIP 4347, HIP 29740, 
HIP 34795, HIP 55852, HIP 58956, HIP 76605 and BD 
+3°2688) have been monitored for many years using the 
CORAVEL spcctrovelocimeter ( Baranne et al.L ri979|) . as 
part of a larger program aiming at finding the frequency 
of spectros copic binaries among s-process-rich late-type 
giant s (see IJorissen fc MavoAll988lll992Hjorissen et all 
Il998l for details and other results from this CORAVEL 
monitoring). Individual radial- velocity measurements for 
th ose stars (in the C ORAVEL-ELODIE syste m as defined 
in lUdrv et al.lll999(l are given in Tables[7|or in lUdrv et al.l 
l|l998(l . For a few other stars (HIP 43042, CS 22942-019 
and CS 22948-027), radial velocities were monitored using 
other inst ruments, and their results were taken frorn the 
literature (|PrestoTi SnedenL l200lt lz"a,cs et a,].[ l200,5|>. 

Orbits were already avai lable for HIP 4347 and 
HIP 29740 JUdrv et all. [l99i), HIP 43042 (IZacs et all 
l2005l). as wen as for CS 22942-019 and CS 22948-027 
(IPreston fc Snedenl. l200ll) . A new orbital solution has 
been derived for BD +04°2466 (Tableland Yig.^. The 
binary nature of those stars is therefore beyond doubt. 



The radial-velocity standard deviation of HIP 58596 
is larger than expected based on the uncertainty on one 
measurement (Fig. EJ, but no satisfactory orbital solu- 
tion could be found. A 3-d orbit (with an eccentricity 
of 0.30) is possible, but this short orbital period is not 
consistent with the giant nature of HIP 58596, which 
imposes orbita l perio ds of at least 20 d (see Fig. 4 of 
IPoiirbaix et al.L l2004ll . The large standard deviation ex- 
hibited by HIP 58596 is therefore very likely another 
example of the large intrinsic jitter of ten obser v ed for 
met al-deficient stars, a s discussed by iMcClurel (|l984l) 
and lCarnev et al ] J2003|) . Nevertheless, the proper motion 
analysis presented in Sect. 14.31 hints at a possible long or- 
bital period for HIP 58596. 

A 14 y radial-velocity monitoring for HIP 34795 with 
the northern and southern CORAVELs is not very con- 
clusive either, mostly because there are difficulties in find- 
ing the zero-point offset bet ween the two instru ments for 
such large radial velocities ljUdrv et al.l Il999j) . When a 
— 1 km s^^ offset is applied to the northern velocities (with 
respect to the values listed in Table |7|), a long-term trend 
seems to be present, albeit with some superimposed jit- 
ter (Fig. [TT|l . The analysis of the Hipparcos astrometric 
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Table 4. Abundances for HD 139409, in the scale where 
loge(H) = 1 2. and normalized with r espect to the solar 
abundances ijGrevesse fc Sauval . Il998|) . The standard de- 
viations a and the number N of lines used in the analysis 
are also given. 



Element (X) 


Z 


loge(X) 


a 


iV 


[X/Fe] 


Cl 


6 


8.46 


0.10 


4 


+0.36 


Nal 


11 


6.02 


0.05 


5 


-1-0.11 


Mgl 


12 


7.31 


0.12 


6 


+0.15 


Si I 


14 


7.33 


0.06 


20 


+0.20 


Cal 


20 


6.20 


0.07 


6 


+0.26 


Sen 


21 


2.79 


0.06 


10 


+0.04 


Til 


22 


4.77 


0.09 


51 


+0.17 


Till 


22 


4.72 


0.11 


4 


+0.12 


Vl 


23 


3.76 


0.11 


23 


+0.18 


Cri 


24 


5.21 


0.11 


11 


-0.04 


Crll 


24 


5.36 


0.09 


3 


+0.11 


Mm 


25 


5.13 


0.11 


9 


+0.16 


Fei 


26 


7.08 


0.08 


136 




Fell 


26 


7.08 


0.12 


14 




Nil 


28 


5.86 


0.07 


44 


+0.03 


Yi 


39 


2.22 


0.10 


2 


+0.40 


Yii 


39 


2.41 


0.07 


6 


+0.59 


Zri 


40 


2.58 


0.09 


3 


+0.40 


Zrll 


40 


2.78 


0.12 


3 


+0.60 


Ban 


56 


1.89 


0.12 


3 


+0.18 


La II 


57 


1.05 


0.12 


3 


+0.30 


Ceil 


58 


1.46 


0.14 


5 


+0.30 


Ndii 


60 


1.38 


0.09 


6 


+0.30 




Table 6. Orbital elements for BD +04° 2466. 



4.6x10-' 4.8x10" 5x101 

(HJD - 2 400 000) 

Fig. 8. Radial velocities as a function of heliocentric 
Julian Day for BD +4° 2466, superimposed on the orbital 
solution corresponding to the orbital elements listed in 
Table El 

The basic idea is to quantify the likelihood of the fit 
of the Hipparcos as trometric data with a n orb ital model. 
For that purp ose, 'Po urbaix fc Arenoul l)200l|) (see also 
Ijancart et all i2005,') introduced several statistical indica- 



P(d) 


4592.7+51.1 


e 


0.286+0.02 


uj (deg) 


266.8+5.3 


Vj (km s^^) 


+39.29+0.11 


K (km s~^) 


5.67+0.12 


T (JD) 


2 445 076.6+73.0 


/(M) (Mo) 


0.076 


oi sin i (Gm) 


343.18 


iV 


41 



data presented in Sect. 14.21 suggests that the star might 
be binary, although the evidence is not very conclusive. 

Finally, there is no sign of radial- velocity variations for 
BD +3°2688 (Fig.UHI). 

4.2. Orbital motion from Hipparcos astrometric data 

A tailored reprocessing of the Hipparcos Intermediate 
Astro metric Data (hereafter IAD; Ivan Leeuwen fc Evansl 
Il998|) makes it possible to look for a possible orbital signa- 
ture in t he astrometric motion , follow i ng the method out- 
lined b y IPourbaix fc JorissenI l)2000() . IPourbaix fc BofSnl 
ll200.'^.IPourbaixl ^l2004^ and applied to barium stars by 
I.Torissen et al.l l(2004lJ) . We give here only a brief summary 
of the method. 
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^ 26t) 
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K 264 
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4.8x10^ 



J , I I !_ 

4.9x10" 
(HJD - 2 400 



5x10" 



000) 



Fig. 9. Radial velocities as a function of heliocentric 
Juhan Day for HIP 58596 = HD 104340 (open circles 
correspond to measurements obtained with the southern 
CORAVEL, and filled circles with the northern one). 
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4.95x10" 
(HJD - 2 400 000) 



5x10* 



Fig. 10. Radial velocities as a function of heliocentric 
Julian Day for BD +3°2688. 



tors which allow us to decide whether to keep or to discard 
an orbital solution. Those indicators relevant to our pur- 
pose are the following: 




312 



4.6x10-' 4.7x10" 4.8x10" 4.9x10" 
(HJD - 2 400 000) 



5x10" 



Fig. 11. Radial velocities as a function of heliocentric 
Juhan Day for HIP 34795 = HD 55496. The symbols are 
as in Fig. |51 With respect to the data listed in Table [TJ 
an offset of —1 km s^^ has been applied to the northern 
measurements to make them consistent with the southern 
ones. 



The addition of 4 supplementary parameters (the four 
Thiele-Innes orbital constants) describing the orbital 
motion should result in a statistically significant de- 
crease of the for the fit of the N IAD with an orbital 
model with 9 free parameters (xt)j e^s compared to a 
fit with a single-star solution with 5 free parameters 
iXs)- This criterion is expressed by an i^-test: 



Pr2 = Pr[F > F(4,7V-9)], 
where 



F = 



N- 



9 Xl 



(1) 



(2) 



Pr2 is thus the first kind risk associated with the rejec- 
tion of the null hypothesis: "t/iere is no orbital wobble 
present in the datd'' . 

Getting a substantial reduction of the with the 
Thiele-Innes model does not necessarily imply that the 
four Thiele-Innes constants A,B,F,G are significantly 
different from 0. The first kind risk associated with the 
rejection of the null hypothesis "t/ie orbital semi-major 
axis is equal to zero" may be expressed as 



Pra = Pr[x 
where 



2 

ABFG 



> x'm, 



Xabfg — X^V ^X, 



(3) 



(4) 



and X is the vector of components A,B,F,G and V 
is its covariancc matrix.^ 

For the orbital solution to be a significant one, its pa- 
rameters should not be strongly correlated with the 
other astrometric parameters (e.g., the proper mo- 
tion). In other words, the covariance matrix of the as- 
trometric solution should be dominated by its diagonal 
terms, as mea sured by the effi ciency e of the matrix be- 
ing close to 1 l)Eichhornlll989() . The efficiency is simply 
expressed by 



n 



(5) 



where Xk and Vkk are respectively the eigenvalues and 
the diagonal terms of the covariance matrix V. 

With the above notations, the requirements for a star 
to qualify as a binary is then 



= {Pr2 + Pr^)/e < 0.02, 



(6) 



where the threshold v alue of 0.02 has been ch osen to min- 
imize false detections ((.Torissen et al.LEo04bl) . 

^ Since it may be shown that xs ~ Xt = Xabfg^ the Pr2 
and Pra tests are in fact equivalent provided that Xt ^ N — 9. 
Thus, if Pr2 and Pra are significantly different, it means either 
that the Thiele-Innes orbital model does not fit the data very 
well (xt >> N ~ 9), or that it fits much better than could be 
expected (xt << A*' — 9). We are indebted to L . Lindegren for 
this clarification (see also Ijancart et alll2005^ . 
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Hipparcos data are, however, seldom precise enough to 
derive the orbital elements from scratch. Therefore, when 
a spectroscopic orbit is available beforehand, it is advan- 
tageous to import e, P, T from the spectroscopic orbit and 
to derive the remaining astrometric elements (as don e by 
IPourbaix fc JorissenL I2OO0I: IPourbaix fc Boffinl . l2003|) . If a 
spectroscopic orbit is not available, trial (e,P, T) triplets 
scanning a regular grid (with 10 < P(d) < 5000 imposed 
by the Hipparcos scanning law and the mission duration) 
may be used. The quality factor a is then computed for 
each trial (e, P, T) triplet, and if there exist triplets yield- 
ing a < 0.02, the star is flagged as a binary. To test its suc- 
cess rat e, this method has been applied bv I Jorissen et all 
ll2004bt) on a sample of barium stars. These authors show 
that, when zu > 5 mas and 100 < P(d) < 4000, the (as- 
trometric) binary detection rate is close to 100%, i.e., the 
astrometric met hod recovers a l l know n spectroscopic bi- 
naries (see also Ijancart et al.l l2005|) . When the orbit is 
not known beforehand, the method makes it even possible 
to find a good estimate for the orbital period, provided, 
however, that the true period is not an integer fraction, or 
a multiple, of one year. Here the method is applied to the 
sample of metal-deficient barium stars listed in Table |5l 

The method flags as definite binaries the stars 
HIP 29740, 34795, 43042, 97874 and 107478 (Figs. [Eland 
[T^ . In two cases (HIP 29740 and 43042), the IAD method 
thus confirms the conclusion from the radial- velocity mon- 
itoring, but yields as well three new binaries (HIP 34795, 
97874 and 107478). Two spectroscopic binaries (HIP 4347 
and HIP 55852) are not detected by the IAD method be- 
cause of their small parallax or long orbital period. The 
non-binary nature of HIP 58596, already suspected from 
the radial-velocity data, is confirmed by the analysis of 
the IAD (Fig.CSI). 

4.3. Orbital motion from a comparison of Hipparcos 
and Tycho-2 proper motions 

an fc Makarovl l(2003() suggested that the comparison 
of Hipparcos and Tycho-2 ( Heig et all l2000b() proper mo- 
tions offers a way to detect binaries with long periods (typ- 
ically from 2000 to 4000 d). The Hipparcos proper motion, 
being based on observations spanning only 3 y, may be al- 
tered by the orbital motion, especially for systems with 
periods in the range of 2000 to 4000 d whose orbital mo- 
tion was not recognized by Hipparcos. On the other hand, 
this effect should average out in the Tycho-2 proper mo- 
tion, which is derived from observations coverin g a much 
longer time span This m eth od, already used bvlMakarovl 
l|2004j) . IPourbaixl l|2004l) and l.Iancart et all !{200^ . works 
best when applied to stars with parallaxes in excess of 
about 5 mas. 

The method evaluates the quantity 



Xobs = (MHIP - MTyc)* W ^ (/XHIP - MTyc), 



(7) 



where /xhip and /xxyc are the vectors of a and S com- 
ponents of the Hipparcos and Tycho-2 proper motions. 



respectively, and W is the associated 2x2 variance- 
covariance matrix. The covariance between fia,iiip and 
/^AHTPi as p r ovided by field H28 of the Hipparcos cata- 
logue |)ESAI . Il997tl and the correlation between Tycho- 
2 and Hipparcos proper n iotions, as encapsu lated in the 
quantity R of Table 1 of lH0g et all l)2000a|) . have both 
been been considered (see Jancart et al. 2006 for details). 

Since the above quantity follows a probability dis- 
tribution function with 2 degrees of freedom, it is then 
possible to compute the probability Prob(x^ > Xobs)> 
giving the first kind risk of rejecting the null hypothesis 
MTycho = /^Hip while it is actually true. This probability 
is listed in Table[Sl along with Xobs^ ^^^"^ ^^'^ fiagged 
as binary if Prob< 0.1. 

Only HIP 29740, HIP 76605 and HIP 97874 satisfy the 
test at the 10% threshold. Note, however, that all the other 
stars have parallaxes smaller than 5 mas, which make the 
test less efficient. 

5. Summary of the binary criteria and discussion 

The situation may be summarized as follows (see also last 
column of Table (S)) : 

— Definite binaries with known orbits: HIP 4347, 
HIP 29740 (passes all three binarity tests), HIP 43042, 
HIP 55852, CS 22942-019, CS 22948-027; 

— Suspected binaries from astrometric data (either IAD 
or proper motions; no or inconclusive radial- velocity 
data): HIP 34795, HIP 76605, HIP 97874 (both astro- 
metric tests yield positive results), HIP 107478; 

— Data inconclusive (mainly because of too small a par- 
allax): HIP 11595, HIP 25161, HIP 69834; 

— Non-binary stars (mainly from radial-velocity 
data, not contradicted by astromctry): HIP 58596, 
BD -h3°2688. 

The latter two non-binary stars are in fact good can- 
didate thermally-pulsing AGB stars, as revealed by their 
location in the HR diagram (Table |2] and right panel of 
Fig. Hence, they must not be binaries. 

Finally, we come to the central question of this paper: 
Why do the metal-deficient barium stars, despite being 
binaries and occupying the same location of the HR dia- 
gram as YSyS, do not exhibit symbiotic activity? Three 
possible answers have been suggested in this paper: 

(i) Some among the stars listed in TableEland displayed in 
Fig. n are in fact not barium stars (especially HIP 80843 
= HD 148897). 

(ii) Among those which arc barium stars, some seem to lie 
on the TP- AGB, and thus need not be binaries (HIP 58596 
= HD 104340, BD +3°2688). They therefore cannot ex- 
hibit symbiotic activity. 

(iii) Finally, there remain a few genuine metal-deficient 
barium stars in the sample. Why are they not symbi- 
otic stars? It seems that the answer to that question lies 
in the different period distributions for YSyS and metal- 
deficient barium stars: YSyS have shorter orbital period 
than metal-deficient barium stars, as seen in Fig. 1141 
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10 100 1000 10 100 1000 
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Fig. 12. The a statistics (Eq. 0) as a function of the trial orbital period (assuming e = 0) for metal-deficient barium 
stars. The period from the spectroscopic orbit (represented by an arrow) of HIP 29740 (upper left panel) indeed lies 
within the range of minimum a values. By comparison, the suspected binaries HIP 97874 and HIP 107478 are likely 
to have periods P ^ 800 d and ~ 250 — 450 d, respectively. For comparison, stars non-flagged as binaries have no 
(e,P, r) grid points with a < 0.02 (see Fig. \Vi\ . The vertical dashed lines represent multiple, or integer fractions, of 
1 y. At those periods, there is a strong correlation between the parallactic and orbital signals, which degrades the a 
statistics and makes binaries difficult to find at those 1-y alias periods. 



This argument seems to apply especially to HIP 29740 
(=HD 43389), which has been assigned a very bright ab- 
solute visual magnitude of —3.5 by the maximum likeli- 
hood method of .Mennessier et al It is therefore 
expected to have a rather strong mass loss rate, and be 
a good candidate YSyS. However, with its orbital period 
of 1689 d, it lies at the long-period edge of the period 



distribution of YSyS (Fig. 114(1 . The same difference seems 
to exist between the pe riod distributions of red s vmbi- 
otics and binary S stars l)Van Eck fc .TorissenL l200^ and 
Fig.lHI). 

However, a firm conclusion on this issue should await 
the determination of the orbital periods for the metal- 
deficient barium stars flagged as binaries by the IAD 
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Fig. 13. Same as Fig.inifor HIP 34795 (left panel, flagged as binary) and HIP 58596 (right panel, non-binary) . 



method (especially HIP 11595, HIP 25161, HIP 69834, 
HIP 97874, HIP 107478), so as to make the comparison 
between the orbital period distributions of metal-deficient 
barium stars and YSyS more meaningful. 
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Table 7. Individual CORAVEL radial velocities and 
associated errors for HIP 34795, HIP 58596, HIP 76605, 
BD +3°2688 and BD +04° 2466. The last column indicates 
which one of the two CORAVEL spectrovelocimeters has 
been used (NO = CORAVEL north at the Observatoire 
de Haute Provence; SO = CORAVEL south on the 
Danish 1.54-m telescope at ESO) 



HIP 34795 = 


HD 55496 






DDMMYY 


HJD 


R.V. 


e 








(km s ) 


/I — 1 \ 
(km s ) 




050183 


45340.695 


315.30 


0.33 


SO 


151283 


45684.797 


314.79 


0.33 


SO 


120684 


45864.434 


315.71 


0.46 


SO 


310185 


46097.658 


314.78 


0.34 


SO 


271185 


46397.826 


314.95 


0.30 


SO 


111186 


46746.858 


314.56 


0.31 


SO 


101286 


46775.595 


316.07 


0.56 


NO 


150187 


46811.682 


314.98 


0.31 


SO 


160288 


47208.394 


316.01 


0.49 


NO 


231088 


47458.704 


315.83 


0.70 


NO 


151288 


47511.857 


315.11 


0.33 


SO 


090289 


47567.393 


315.39 


0.56 


NO 


241189 


47855.642 


315.44 


0.57 


NO 


070390 


47958.339 


317.36 


0.66 


NO 


050490 


47987.497 


314.78 


0.31 


SO 


061190 


48202.697 


315.63 


0.45 


NO 


031290 


48229.645 


316.10 


0.56 


NO 


180391 


48334.380 


314.95 


0.62 


NO 


110393 


49058.361 


315.00 


0.67 


NO 


130393 


49060.338 


314.32 


0.80 


NO 


180194 


49371.471 


315.37 


0.58 


NO 


301295 


50082.752 


313.99 


0.34 


SO 


231196 


50411.600 


316.66 


0.56 


NO 


120197 


50461.480 


317.72 


0.63 


NO 


HIP 58596 = 


HD 104340 






DDMMYY 


HJD 


R.V. 


e 








(km s '-) 


(km s ^) 




141289 


47875.832 


265.29 


0.45 


SO 


140392 


48696.805 


264.18 


0.30 


SO 


120393 


49059.755 


266.40 


0.29 


SO 


190393 


49066.707 


264.89 


0.33 


SO 


050493 


49083.708 


265.41 


0.28 


SO 


200493 


49098.453 


263.93 


0.46 


NO 


080394 


49420.525 


263.34 


0.47 


NO 


030494 


49446.435 


266.38 


0.41 


NO 


140594 


49487.360 


265.37 


0.46 


NO 


260195 


49744.638 


264.59 


0.56 


NO 


300395 


49807.475 


263.16 


0.52 


NO 


210795 


49920.516 


264.34 


0.30 


SO 


191295 


50071.721 


265.60 


0.37 


NO 


150296 


50129.556 


263.70 


0.43 


NO 


130197 


50462.684 


265.61 


0.45 


NO 


HIP 76605 = 


HD 139409 






DDMMYY 


HJD 


R.V. 


e 








(km s^^) 


(km s^^) 




090790 


48082.629 


62.43 


0.30 


SO 


260794 


49560.513 


63.36 


0.27 


SO 
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Table 7. Continued. 



HIP 


55852 = : 


BD +04°2466 




DDMMYY 


HJD 


R.V. 

(km s'^) 


e 

(km s^^) 




070282 


45008.814 


38.11 


0.69 


SO 


130684 


45865.491 


44.73 


0.74 


SO 


030285 


46100.832 


45.14 


0.68 


SO 


290385 


46154.675 


44.78 


0.58 


SO 


200186 


46451.831 


43.54 


0.53 


SO 


140286 


46476.783 


43.32 


0.49 


SO 


291186 


46764.706 


42.05 


0.96 


NO 


090187 


46805.864 


41.36 


0.49 


SO 


080287 


46835.808 


41.99 


0.41 


SO 


280188 


47189.812 


40.50 


0.43 


SO 


040388 


47225.726 


39.83 


0.41 


SO 


250488 


47277.657 


40.17 


0.48 


SO 


051288 


47501.707 


37.84 


0.69 


NO 


060389 


47592.565 


38.14 


0.64 


NO 


220589 


47669.526 


37.64 


0.49 


SO 


241289 


47885.857 


37.48 


0.45 


SO 


050290 


47928.800 


36.83 


0.48 


SO 


160290 


47939.743 


37.05 


0.43 


SO 


260290 


47949.595 


37.00 


0.69 


NO 


040490 


47986.631 


37.01 


0.42 


SO 


110490 


47993.594 


35.99 


0.49 


SO 


050590 


48017.647 


36.35 


0.61 


SO 


020690 


48045.501 


36.82 


0.46 


SO 


260291 


48314.487 


34.56 


0.73 


NO 


260391 


48342.724 


34.92 


0.47 


SO 


150392 


48697.747 


34.26 


0.42 


SO 


120393 


49059.726 


33.77 


0.39 


SO 


070493 


49085.670 


33.15 


0.40 


SO 


080493 


49086.448 


32.60 


0.70 


NO 


100294 


49394.609 


34.54 


0.77 


NO 


020394 


49414.713 


35.84 


0.43 


SO 


070394 


49419.515 


34.92 


0.52 


NO 


080494 


49451.424 


37.51 


0.88 


NO 


140594 


49487.351 


36.78 


0.64 


NO 


310195 


49749.592 


39.24 


0.74 


NO 


150395 


49792.681 


41.24 


0.44 


SO 


020495 


49810.460 


40.55 


0.44 


NO 


190795 


49918.506 


41.94 


0.53 


SO 


191295 


50071.715 


43.32 


0.60 


NO 


130197 


50462.674 


44.33 


0.90 


NO 


200199 


51199.634 


41.42 


0.80 


NO 


BD +03°2688 


DDMMYY 


HJD 


R.V. 

(km s^^) 


e 

(km s^^) 




080493 


49086.481 


33.16 


0.43 


NO 


050394 


49417.506 


32.35 


0.47 


NO 


030494 


49446.510 


32.33 


0.44 


NO 


140594 


49487.388 


33.00 


0.47 


NO 


030295 


49752.658 


33.52 


0.52 


NO 


300395 


49807.506 


32.64 


0.71 


NO 


201295 


50072.729 


33.07 


0.46 


NO 



